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SUMMARY 

The optimal control of large-angle rapid maneuvers and vibrations of a 
Shuttle-mast-reflector system is considered. The nonlinear equations of motion are 
formulated by using Lagrange s formula, with the mast modeled as a continuous beam. 
The nonlinear terms in the equations come from the coupling between the angular 
velocities, the modal coordinates, and the modal rates. Pontryagin’s Maximum 
Principle is applied to the slewing problem, to derive the necessary conditions for 
the optimal controls, which are bounded by given saturation levels. The resulting 
two-point boundary- value problem (TPBVP) is then solved by using the 
quasilinearization algorithm and the method of particular solutions. In the 
numerical simulations, the structural parameters and the control limits from the 
Spacecraft Control Laboratory Experiments (SCOLE) are used. 

In the two-dimensional (2“D) case, only the motion in the plane of an Earth 
orbit or the single-axis slewing motion will be discussed. The effects of the 
structural offset connection, the so-called axial shortening, and the gravitational 
torque on the slewing of a flexible spacecraft in Earth orbit is considered. In the 
case of three-dimensional (3— D) slewing, the mast is modeled as a continuous beam 
subject to 3-D deformations. The numerical results for both the linearized system 
and the nonlinear system are presented to compare the differences in their time 
responses . 


INTRODUCTION 

Future space missions require large-angle rotational (attitude), 3-dimensional 
maneuvering ("slew”) of a large flexible spacecraft in target acquisition, target 
tracking, and surveying multiple targets, etc. The whole spacecraft system may be 
composed of multibodies, including the rigid 1 parts involving large relative 
movements and its flexible parts undergoing ’’small” deformations. The motions of the 
system for these space activities are best described by nonlinear equations instead 
of the linearized” or linear equations. Many authors have considered the problem of 
large-angle rapid maneuvers of flexible spacecraft (refs. 1-8). The direct 
application of Pontryagin’s maximum principle to this problem has been applied by 
many authors (refs. 1-3, 6-8). Most of these researches are concentrated on the 2-D 
slewing problem, except ref. 8, a recent result for the 3-D slewing of the SCOLE. 


2-D Slewings 

In ref. 6, the rapid slewing of the 2-D SCOLE has been considered. It is 
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observed that (ref. 6) the time response history of the nonlinear system has a shift 
from that of the linearized system; the reason for this is mainly due to the 
structural offset (ref. 7). 

The so called axial shortening effect of a beam induced by its transverse 
displacement has been brought to attention by some authors (Refs. l - 3, 9). Although 
the shortening terms have been included in the equations (Refs. 1“3) , their effect 
on the slew lacked quantitative analysis; specifically, the numerical examples with 
and without these terms were not provided. On the other hand, a numerical example in 
Ref. 9 shows that large differences do result between models with and without the 
shortening effect. But the numerical example is only for an uncontrolled dynamical 
response case and the main body’s motion is prescribed. In ref . 7 , therefore, the 
shortening terms are considered in the formulation of the equations of motion and 
numerical examples both with and without these terms are presented to compare the 
difference between them. Also in ref. 7, the gravitational torque terms are modeled 
and included in the equations to show their effect on the slewing motion. 


3-D Slewings 

The direct solution of the open-loop TPBVP for 3“D slews of flexible spacecraft 
resulted in numerical problems with rank— deficient matrices as stated in ref. 5. 
However, a different numerical method may be used to overcome this difficulty. In 
ref. 8, the problem has been solved successfully by using the quasilinearization 
algorithm and the method of particular solutions for 3-D slews of the asymmetrical 
flexible SCOLE configuration. 

The open - loop slewing approach has several obvious distinct properties. First, 
the control law is easy to implement in practice for both ground tests and space 
flight tests. Second, the open-loop solution may serve as a good reference for the 
feedback control law design, as proposed in refs. 4-5, in which the open-loop 
solution for a rigid (instead of a flexible) spacecraft is used as the nominal 
reference trajectory. As an extension to refs. 4—5, it may be helpful if the 
open-loop solution for the 3-D slew of a flexible spacecraft system could also be 
used as a nominal reference solution. In addition, through the present study, we can 
also see how different are the responses of the nonlinear system from those of the 
linearized system. 

In the present report, we will summarize most of the results obtained in refs. 
7-8. At the same time, the detailed numerical techniques (briefly mentioned in refs. 
7 — 8 ) for solving the nonlinear TPBVP will be discussed. Numerical examples are 
presented to illustrate the use of the techniques and the numerical problems 
associated with the calculations will be discussed. 


EQUATIONS OF THE SYSTEMS 
System Configurations 

The 2-D and 3-D models of the orbiting SCOLE are shown in figures 1 and 2 , 
respectively. The Shuttle and the reflector are assumed to be rigid bodies. One end 
of the flexible mast is fixed to the Shuttle at its mass center, o g , while the other 

end is firmly connected to the reflector at an offset point, a f (x f in the 2-D 
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case). Three Euler angles (0^, B^ f 0^) (0 in the 2-D case) or four quaternions are 

used to describe the attitude of the Shuttle with respect to an orbital reference 
system. 

The 3-D deformation of the mast consists of two bending deflections l/(z,t) and 
V(z,t) in the x-z and y-z planes, respectively, and torsion <J>(z,t) about the z axis 
It is assumed that these deformations are small as compared with the length of the 
mast and can be expressed by the following modal superposition formulas (ref. 10): 


U(z,t)=^ i (z)a. (t) , V(z,t)=p 7 j (zja^t), <J>(z,t)=^ i (z)« 4 (t) , (1) 

where , rj^ , and £, are modal shape function vector components normalized by a 
common factor, and is a scaled modal amplitude associated with the ith mode. In 
the 2-D case, only the first equation is used. 

The free vibration of this structure can be considered as a space free-free 
beam vibration problem with boundary conditions including the masses and moments of 
inertia of the Shuttle and the reflector. The partial differential equation 
formulation for this problem (refs. 10-11) can be solved by using the separation of 
variables method. Note that the natural frequencies and modal functions of the 2-D 
structure are different from those of the 3-D structure. 

2-D Dynamical Equations 

After developing the kinetic energy and potential energy, we can obtain the 
dynamical equations of the 2-D structure in the following matrix form (ref. 7): 


I+2a T m a \ 

(m +M a) T 

3 2 : 

2 4 

m +M a 

M 

2 A 

3 


-26a 7 (m +M ct)-« T M a ~6V JdO +Q_ 

3 2 5 y t/ 

-2dM a+£ 2 (m +M a)-Kor-0V /da +Q 
3 a 2 a 


( 2 ) 


where a is the nxl modal amplitude vector, n is the number of flexible modes; 

M. (i=2,3,4,5) and M are constant nxn matrices; m and m are nxl constant vectors; K 
i a a 2 

is the nxn constant diagonal stiffness matrix; V is the gravitational energy of the 

9 

system which is a function of the orbital rate, <i> o , rotation angle, 0, and the 
deformation amplitude, a(t); and Q^, are the generalized forces produced by the 
controls associated with 0 and a, respectively. 

In equation (2), the elements of the vector, m , and the matrices, M , M , and 

a a A 

M have a common factor x , the offset. In another words, x =0 results in m =0, etc. 
5 r r a 

The effect of this structural offset on the slewing is analyzed by changing the 

value of x (from 0 to 32.5 ft). Of course, different x imply different natural 
r r 

frequencies and modes. The frequencies decrease as the offset distance increase. 

The axial shortening of the beam due to the geometric deformation is also 
considered by adding a shortening term 
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Lz- 


4r $ 


dz 


into the formulation. Apparently, Az is a second order term in the modal amplitude, 
a. The shortening term is involved in the matrices , M 5 in equation (2). 

If we distinguish the terms on the right side of equation (2) according to 
their order in a and a, we can see that m g and M a are involved in the lower order 

terms, while M and M contribute to the higher order terms. Therefore, for moderate 
2 5 

nonzero values of x and small-deformation slews, the influence of the structural 

offset can be greater than the shortening effect. The linearized equations can be 
obtained by neglecting all nonlinear terms, 


T : T 
1 i m 2 
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+Q e 

,v 3 
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K 
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? 

i 

+ V 


(3) 


LIN 


where "LIN" refers to constant and linear terms. Note that on the right side of 
equation (3), both the structural offset and the shortening terms disappear. This 
means that by using the linearized equations of motion, we may lose some important 
information such as the terms associated with the structural offset about the 
original system. 


3-D System Equations 

The dynamic equations for the 3-D system can be obtained, by using the 
Lagrangian method, in the following state form (ref. 8); 


0) 

P 


= (A +B )u +(C 0 )u +Doc +(E +F )u 
a p ck 


(4) 


where a is the nxl modal amplitude vector, is the modal rate vector, u> is the 

3x1 angular velocity vector of the Shuttle, w is a 6x1 vector defined by 


o, 2 o ) 3 ^ mwj 


and u is the 9x1 control vector: 

u = [ f f f if f If f If f ]! 

1 lx ly 1 z 2x 2y ' 3x 3y ' 4x 4y 

A , D , and E , are constant matrices. Also, (B ) 

(3+n)x6' ( 3+n ) xn (3+n)x9* « <3+n)x6 

, are defined as 

v 1 - yiv> ! v 1 c 3 pl> V tF ,“ ! v 


(n ) , (F ) , are defined as 

V p ( 3 + n ) x3 Ot^O + nJxS 


B =[B a I B a 
a l 2 


p ,«] 
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where B_ , C, , F, are all (3+n)xn constant matrices. 

Note that in the development of equations (4), a constant inertia (mass) matrix 
has been assumed for convenience. Meanwhile, all the second and higher order terms 
of the flexible variables (a and a) are abandoned in the final equations. But all 
nonlinear terms (o) representing the rigid body motion are retained. 

Clearly, the dynamic equations for the rigidized spacecraft can be obtained by 
deleting all terms related with a and 0, that is, 

G) = Ag) + Eu (5) 

where A and E are 3x6 and 3x9 constant matrices, respectively. 

A linearized form of Equation (4) can also be obtained by deleting all 
nonlinear terms. 


y = Da + Eu ( 6 ) 

By using the quaternion vector q=[q Q q ^ q z qj 1 , the kinematic equations can be 
expressed as 
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OPTIMAL CONTROL-NECESSARY CONDITIONS 

For the 3-D slewing problem, we use the following quadratic cost functional, 

IP 1 

J= 5J f ^ CfTQ 1 a '*w T Q 2 » +P T %P +u T Ru)dt (8) 

where , and R are weighting matr 
of the controls, u, are bounded, 

hi* 

The Hamiltonian of the system 
H ~(a T Q a +o T Q (o +0 T Q 0 +u T Ru 

where p, y, and \=[\ i A 2 ] T are the 
respectively. By using the Maximum 

unrestricted optimal control problem are the dynamical equations (4, 7) plus the 
following differential equations for the costates. 


i 


ices, t f is the given slewing time. The magnitudes 


i b 

is, 


, 9. 


(9) 


+p uq)+y T 0+X T [(A +B^)(i) +(C ft )w +Da +(E +Fju] 


P' 


a 


costate vectors associated with q, a, u, and 0, 
Principle, the necessarv r.ondi tinnc 
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(10a, b) 


p =-t? = 1 « p. i =-§5 = “V -d t x -(£x)5 -(f>)u 


8 q 2 — 
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i.Hs? = -Q^tqJp -[x T (A + i a )](o -(c^) T x, x 2 =-|g = -y -(^x)« (io c ,d) 


1 6 co "22 
and the optimal control, 


§5=0, ~> U=-R-<E.P/X 


(ID 


The control law (11) is then modified by consideration of the saturations (9). 

otherwise , 
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( 12 ) 


The same type of cost functional has been used for the 2-D slew problem. Therefore, 
the associated necessary conditions are similar to equations (10) and (12). 


Two-Point Boundary- Value Problem (TPBVP) 


One way of obtaining the optimal control law is to transform the above 
necessary conditions into the following TPBVP. Let x represent the state vector, and 
X represent the costate vector. After substituting the control expressions (12) into 
equations (4) and (10b), one can obtain two sets of ordinary differential equations 
for the states and the costates , 


i= V X ’ X) (7*2n)*l 


X f 2 (X ’ X \? + 2n)xl 


(13a) 

(13b) 


with the following boundary conditions, 

x(0) and x(t^) prescribed, X(0) and X(t f ) unknown. 


(14) 


Due to the known boundary conditions being specified at the two ends of the slewing 
time, this problem is usually called the two-point boundary-value problem. This kind 
of split boundary conditions usually result from the large-angle maneuver 
requirements, in which the initial (t = 0) and final (t-t^) states of the system are 

specified. By solving this problem, we can obtain the optimal control (based on the 
necessary conditions). The often used solution strategy is to change the boundary 
value problem to the initial value problem, i.e., find \(0), the missing initial 
costates. Once \(0) is obtained, one can solve the equations (13) as an initial 
value problem by using any numerical integration method. However, owing to the 
nonlinearity of the equations, there is generally no analytical solution to this 
problem or simple numerical method to obtain the solution except for some very 
simple cases such as the linear time— invariant case. The numerical iteration method 
is the general approach to the this problem. 

To start an iteration process, one usually needs an initial guess of X (0> (0). 
Then, equations (13) or their equivalent form (the linearized version of (13)) are 
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solved and a x ( *(t ) i s obtained. Based on the difference Ax(t )=x* 0, (t )-x(t ), 
r f f f 

the correction to X(0), AX(O), is obtained. This gives us a new initial value of 
X(0), X (0). Hen^e^ the ne^t^ iteration begins. The iteration process can be 
terminated when |X (0)~X (0)| is less than a given error limit. One can see 

immediately that if the beginning guess X* * (0) is close to the true value 
(converged value) of X(0), the solution will converge and less iterations are 
needed. However, a "good" guess of X(0) is often difficult to obtain for the general 
nonlinear problems. 

Therefore, the effort for solving the TPBVP is two fold. The first is try to 
establish a good iteration (correction) method with a wide convergence interval so 
that it can guarantee convergence even for a "poor" initial guess. The other is try 
to find a "good" initial guess based on the characteristics of the practical problem 
and using some simplified mathematical models. In this report, we use the 
quasilinearization method. We also use the solution of X(0) from the simplified 
linear, time - invariant model of the system as the initial guess for starting the 
iteration process. 


Linear and Time-Invariant TPBVP 

For linear, time-invariant versions of equations (13) (refs. 1-3), 

z=Az , where z T = [ x T , X T ] 

its transition matrix (constant exponential matrix), 


(15) 


At 


A A 

11 12 

A A . 

21 22 J 


can be used to obtain the initial costates (closed form solution): 

X(0)=A^[x(t.)-A_x(0)] 


11 


(16) 


Nonlinear TPBVP 


The continuation (relaxation) process (to increase the participation of the 
nonlinearity in the solution) and the differential correction (for determination of 
the initial costate variables) have been used in references 1-3 for the 2-D slewing 
problem. However, as stated in ref. 5, the extension of these techniques to the 3-D 
slewing problem has encountered a numerical problem: rank deficiency. 

In references 6-8, the quasilinearization method has been successfully used. 

In this method, one needs to linearize the differential equations (13), 

i=g(z), where z T =[x T , X T ], g T =[fJ, f*] (17) 

about an approximate solution of this equation in the following form (a series of 
linearized, time-variant, nonhomogeneous equations): 

z(l, + 1) “(dg/3z)z (k * 1> +h(z lk) ) (18) 


4-00 


. Ck) 

where z 

approximate 

conditions , 

equations , 

12 ): 


is the kth solution of the same linearized equation. It is also the kth 
solution of the original nonlinear equations (17). Here, the boundary 
(14), are naturally adopted as the boundary conditions of the linearized 
(18). The control expressions, (11), also need to be linearized (ref. 


u (k + l) =u (k)_ R -l [ p (Aa)] T x (k)_ R -l [E+ p (a Ck) ) jT AX (19) 


where Aa=a 


(k + l) 


<k) 


nx 


and AX=X (k+1> -X (k> . By assuming that 


u (k, =-R- 1 [E + F(a (k) )] T X (k) 

for the unbounded control case, equation (19) can be rewritten as, 

u (k+ 1 ) =-R- 1 [F(Aa)] T X lk, -R- 1 lE + F(a <k) )] T X <k ' 1) 

However, in the bounded control case, equations (12) are considered, that is, 

or u 


( 20 ) 


( 21 ) 


(k) 


-u 


lb 


ib 


( 22 ) 


-{R- 1 [E + F(a (k, )]V k, } 1 


Accordingly, at the (k+l)st step, u (k+1) can be determined by 


( k + 1 ) _ 
u 

i 


-u or u , if |{R' 1 [E+F(a (k) )]V k) } u 
ib ib 1 1 1D 

-{R' 1 [F(Aa)] T X (k) -R' 1 [E+F(a <k) )] T X <k+1) } 1 


(23) 


So far, an iteration process is formed. In each iteration, only a linear TPBVP is 
solved. It is this property that gives this approach the name quasilinearization 
method . 


The linear TPBVP can be solved by many ready-made methods. One of 
frequently used algorithms is the method of particular solutions (ref. 
represent the number of the states (also the costates). Equations (18) 
rewritten in the following form. 


the 

13). Let m 
can also be 


x(t)=G(t)x(t)+H(t)X(t) , X(t)=I(t)x(t)+J(t)X(t) 


(24a, b) 


Prom the theorem of the linear system, any solution equation (24a) can be expressed 
as the linear combination of its m+1 particular solutions, i.e., 

m+1 nil 

x(t)= E c x‘(t), as long as E c =1 (25) 

i ■ 1 i 1,1 1 

where c are constants and x‘(t) are the ith particular solution vectors. The method 

begins by integrating equations (24a, b) forward mr+T times, with the initial 
conditions , 
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z 1 (0)= 

’ x (0) 
1 
0 
0 

N 

X“N. 

o 

ii 

‘ x(0) 
0 
1 
0 

z”(0)= 

x(0) 

0 

0 

, and z m+1 (0)= 

“ x(0) 
0 
0 
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- 0 - 


• 0 - 


0 

L i -1 


■ 0 - 


This gives us m+1 particular solutions, x 1 (t), x 2 (t), •••, x m+1 (t). Substituting 
these solutions into equations (25), and setting t=t f , we have 


■ ♦1 m+ 1 

i«i c i xl(t r )=x(t r ) ' i5i V 1 


(26) 


This is a set of m+1 algebra equations for m+1 unknown constants, c . By assuming 

i 

the existence of inverse of the coefficient matrix, we can obtain the solution, 
c= f c , c , c ] T and c ..By doing the following manipulation, 

i ^ n m+ 1 


m+1 

z(0)= E 
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L Rl J 


one immediately realizes that c=X(0), the missing initial costates. 


NUMERICAL RESULTS 
2-D Slews 

The parameters of the orbiting SCOLE (refs. 10-11) are used. The orbital 
angular rate is chosen as u o =0.001 (rad/s) (orbital altitude b« 981 km). The first 

three natural frequencies used here are: 

0.3365257, 2.062547, 5.316669 (hz), for x =0; 

0.3199540, 1.287843, 4.800169 (hz), for x -32.5 (ft) (same as ref. 10) 

All simulations are 90 degree rest-to-rest slews and can be represented by: 

Case 1 x r =0, u^ , R=lE-6, t f =27.6 (s) 

Case 2 x =0, G=[u u u u 1 T , R=DIAG(lE-6, .15, .21, IE-4), t =8.196 (s) 

Case 3 X r =32.5 (ft), u=u^ , R=1E~6, t =27.6 (s) 

Ca»e_4 x f =32.5 (ft), G=(u i U;J ] T , R=DI»G(lE-6, .15, .21, IE-4), t =8.196 (s) 


<h02 



Figs. 3a-g display the time histories of 0(0. u(L,t), ♦(L.t). u(t) for Case 4. 
Clearly, the response of 9(0 for both linear and nonlinear systems are very clos . 
However, there exist some differences between the two systems in u L t). ♦(L.t) and 
the controls, u. The difference is primarily due to the offset ( ’ r 

ft). When x -0, this difference can be reduced markedly, regardless of whether the 

shortening effect and gravitation are considered. It is also interesting to know 
that the controls have large differences only around the mid-slew time. 

Table 1 lists the maximum (minimum) values for the displacement, u(L,t), and 
angle <J»(L,t), of the beam during the associated slews for all cases. The firs m 
in each case lists the results for the linearized system, while all remaining lines 
represent those for the nonlinear system with different considerations. For exampl , 
AL=0 means the shortening effect is not considered. The last column gives t e 
largest relative displacement error, with respect to the linear results, based on 


(e 


D i s p 


),=Max[|MAX.-H*X ulI |/|MAX LI> | , I , „ | / | MIK, , , | ) 


Nonlinear Svstem vs. Lin earized System First, let us examine line 1 and line 2 m 
each case. In Case 1, since no offset, the differences between the two lines are 
very small. In Case 2, where more controllers are used and the slewing time is 
shortened, the differences increase symmetrically ( | MAX| -| MIN| ) , in spite of X f . 

Case 3 uses the same slewing conditions as used in Case 1, except x=32.5 ft. This 

offset shifts the envelop of the response downward and results in a larger relative 
displacement error than that in Case 2. Case 4 is the combination of Cases 2and3. 
The shift now is upward which is due to the inclusion of more controllers. When more 
controllers are used (Ref. 3), the phase of the response reverses, so do the maximum 
(minimum) amplitudes . 

Shortening Effect By comparing line 2 and line 3 in each case, we can see that the 
shortening terms (1) reduce the amplitude (Cases 2 and 3); (2) increase the 
amplitude (Case 4); and shift the response upward (Cases 3 and 4). These 
observations coincide with the fact that Az only results in a second order effect as 
compared with the offset effect. 

Gravitational Effect By observing lines 3 and 4 , we can conclude that the addition 
of the gravitational torques into the equations of motion has a very small effect on 
the slew, although they shift the response downward. This is because the orbita 
rate is much smaller than the slewing rate and the magnitude of the gravitational 
torque term is much smaller than that of the active control torque term. 


3-D Slews 


The location of the mass center of the reflector is X r -18.75 ft, and y f 32.5 
ft. The first five natural frequencies are (hz): .2740493, .3229025, .7487723, 
1.244013, 2.051804. 

The numerical tests based on the previously described method have been 
performed for the roll-axis slews, pitch-axis slews, as well as arbitrary axis 
slews. All these tests are rest-to-rest slews and the iteration process is 
terminated after the initial costates are reached within five digit accuracy. 

The following procedure is designed to obtain the solution of the nonlinear 
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TPBVP (figure 4). First, the linear TPBVP based on equation (6) is solved and a 
nominal trajectory is produced, in which the control is unbounded and the initial 
costates are calculated by using the transition matrix method. Then, a converged 
solution for the linear TPBVP with bounded controls is obtained by iterations 
starting from the previously obtained trajectory. Note that the Euler angles are 
used in all the above computations. 


Next, to obtain the starting solution for the nonlinear TPBVP, the 3 Euler 
angles and the 3 associated costates are transformed to the 4 quaternions and their 
4 costates (from t=0 to t=t f ). Reference 14 provides us the following relationship 

between the quaternions, q(t) and their costates, p(t): 



(27) 


where d are constants. For the case q(0)=[l 0 0 0] T , we can choose p (0)=0, d =0. 

oo 

Then , 


tp,(0) P 2 (0) P 3 (0)]=[d i d 2 d 3 l=d T (28) 

The vector d can be determined by 

d- 2 [initial Euler angle costate vector] (29) 

This result can be proved if one compares the related state and costate equations 
for both linear and nonlinear TPBVPs (for the case (o(0)=0). After finding the q(t) 
by using a nonsingular transformation between q(t) and the Euler angles, 0 (t), 

0 2 (t), and 0 3 (t), one can use equations (27-29) to obtain p(t). 


Finally , the nonlinear TPBVP is solved through the quasilinearization process and 
the method of particular solutions. 

Case 3: Fig. 5 shows the results for a simultaneous three-axis slew (0 =60 

l 

0 2 -3°, 0 3 ~45 deg). The weightings for the states are Q^Q^Q^O. In this case. The 
Shuttle torques (f , f , and f ) and the reflector forces (f , and f ) are 

lx * z 4x 4y 

used. The associated weighting for the control is R=DIAG(lE-4, IE-4, IE-4, 0.6, 

1.4 3). The slewing time, t^ , is 40 sec. The solid lines in the figures 5a— h 

responses of the rigidized system, equation (5), while the dotted lines represent 
the slew results for the flexible spacecraft. 


CONCLUSION 

Generally, for the 2-D case, the linearized system can predict the system 
dynamics very well in the slow slewing case. However, in the rapid large-angle 
slewing problem, the responses of the system deviate noticeably from those described 
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by the linearized equations if the effects of structural offset and axial 

ss.si.-5: s s s«-ir,ssisjs!“^ -- 

neglected in the slew motions considered here. 

The application of Pontryagin's Maximum Principle to the large angle slewing 
maneuver problem has been extended to the slewing of a 3-D flexible spacecra^ 
(SCOLE). A numerical simulation procedure based on t e ^uasi t j 

for solving the resulting nonlinear TPBVP has been established and tested 
successfully for several examples. The general nonlinear dynamica equ 
developed here contain all quadratic terms in the an^lar velocity Wjjnt. of tta 
main body and their coupling with the first order modal amplitudes and ™ dal 
It is suggested that higher order terms be included if a further analysis 
conducted The numerical results show an important fact that the linearized system 
can represent the nonlinear system adequately for predicting the major motions but 
not as well for the secondary motions. For further research, it is reco^ended that 
the applicability of this method to more complicated systems be establis 
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Fig. 3. Results for the 2-D Slew, Case 4: x =32. 
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Table 1 Tip Displacement and Tip Angle 
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SOLUTION PROCEDURE 
FOR NONLINEAR TPBVP 


Linear TPBVP 
Euler Angles 
Unbounded Control 



Linear TPBVP 
Euler Angles 
Bounded Control 


L - 


Transformation: 

Euler Angles -•) Quaternions 
(using relation between them) 

Euler Angle Costates ■ ■) Quaternion Costates 
(using relation between 


quaternions and their costates) 



I l 

Nonlinear TPBVP 


Fig. 4. Solution Procedure for Nonlinear TPBVP 
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(i) U(L,t), V(L,t), Bending (ft) 


V V\^ rH 




U, IN X-Z PLANE 

V, IN Y-Z PLANE 


Time (sec) 

(j) Phi(L,t), Z-Torsion (deg) 



Time (sec) 

Results for the 
Shuttle Torques 
R=Diag(lE-4, IE- 

3-D slew. Case 3: 60-30-45 deg-Slew 
+ Reflector Forces, 

4, IE-4, 0.6, 1.4E-3). 






